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I. INTRODUCTION 



In a series of papers [1, 2] the dynamical equations for elastic deformable astronomical 
bodies in the framework of Einstein's theory of gravity (GRT) were derived. Fields of 
application range from problems of seismology of astronomical bodies including the problem 
of normal modes to the global motion of such bodies in space. Actually one of our main 
interests is the problem of global geodynamics where the accuracy of modern geodetic space 
techniques such as VLBI has long reached the level where effects from relativity have to be 
taken into account. Unfortunately the field of relativistic global geodynamics has not yet 
reached a satisfactory level despite the fact that published amplitudes of nutation series are 
of //as (microarcsec) accuracies. 

This paper presents another necessary step to improve this situation. It extends our pre- 
vious work mentioned above that was based upon two frameworks: (1) the Cartcr-Quintana 
formalism [3, 4] for the description of elastic deformable bodies in GRT by means of a dis- 
placement field and (2) the D amour- Soffel-Xu one (the DSX scheme [5, 6, 7, 8]) on relativistic 
celestial mechanics in the first post-Newtonian approximation to GRT. 

Here we shall expand the general relativistic perturbation equations for the dynamics of 
elastic deformable bodies in terms of Generalized Spherical Harmonics (GSH). This turns 
the set of partial differential equations for scalars, vectors and tensors into a set of ordinary 
differential equations. Actually in the Newtonian approximation numerical programs that 
have been written mainly to deal with the global dynamics of the Earth (e.g., Wahr [9], 
Dehant [10], Schastok [11]) by considering its elastic properties in a local framework usually 
integrate such (Newtonian) equations. In general relativity, the general expansion for vectors 
and tensors by means of GSH has been discussed in several papers [12, 13]. Here, we 
apply these methods and expand the general relativistic dynamical equations for elastic 
deformable, nonrotating astronomical bodies by means of GSH. 

In the following we consider the dynamics of one of these bodies in its own local coordinate 
system (cT, X) restricting ourselves to a spherically symmetric non-rotating relaxed ground 
state. The extension to some rotating axially symmetric ground state with nonvanishing 
dynamical ellipticity so that problems of precession and nutation can be treated will be the 
subject of another paper. In Sec. II, we briefly review the main results in Ref. [1]. We also 
rewrite the equations and introduce suitable notation in order to facilitate the comparison 
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with the corresponding Newtonian equations (Ref.[9]). 

In Sec. Ill, Generahzed Spherical Harmonics are introduced. Our main new results are 
presented in Sec. IV, where the post Newtonian perturbation equations (PDE) are expanded 
so that they reduce to a set of ordinary differential equations (ODE). In Sec. V some con- 
clusion can be found. 



II. PERTURBED FIELD EQUATION AND PN EULER EQUATION 

Let us first recall some relations that will be relevant for this paper. In the DSX- 
formalism, the Einstein field equations to first post-Newtonian order can be written as 

[5] 

^^-^^ = -47rGE + 0(4), (2.1) 
V^W^" = -47rG'S" + 0(2) , (2.2) 

where W and W"" are the scalar and vector potentials that describe the gravitational inter- 
action, and the gravitational mass-density and mass- current-density, S and S", are related 
with the energy momentum tensor by S = (T°° -|- T"")/c^, = T^^/c. We shall often 
abbreviate the order symbol 0(c~"') simply by 0{n). 

For a non-rotating static spherically symmetric body (unperturbed state), the energy- 
momentum tensor takes the form 



- pc^ (l + ^) + 0(2) ' (2-3) 
T<^« = 3p(^l-— j+0(4), (2.4) 
rO« = 0(3) , (2.5) 



where pc is the energy density and p is the pressure. The field equations reduce to 



+ 0(4) , (2.6) 



V^iy" = 0(2) , (2.7) 

where p* = p + ^ is the chemical potential per unit volume. For a static equilibrium 
configuration the scalar potential 1^ is a function of the radial coordinate r only, and the 
gravito-magnetic potential W"" vanishes. 
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Now we consider a perturbed state. The perturbed field equations for the Eulerian 
variations of gravitational potentials, SW and SW"' read 

1 d^6W 



VW---^ = -47rG'(5E + 0(4), 
V^SW = -47rG'5E" + 0(2) , 

where, by using Eqs.(4.27), (4.31), (4.40) and (4.41) of our previous paper (Ref.[l]), 

= 6p+\ {2p6W + 2W5p + 35p) 



(2.8) 
(2.9) 



- -V • (ps) - 

51T = ps" + 0(2) 



V • (ps) + Bps • VW + pm + 2V • (piys) + 3kV • s] + 0(402,. 10) 

(2.11) 



and K is the compression modulus. Then the field equations can be expressed as 
V''5W - = 47rG(v • (ps) + ^ [v • (ps) + 3ps • VW + p5W 

+2V • {pWs) + SkV • s] I + 0(4) , 
V'^dW = -47rGps" + 0(2) . 



(2.12) 

(2.13) 



The perturbed PN Euler equation (for non-rotating, static and spherically symmetric ground 
state) reads (Eq.(4.32) of Ref. [1] for Q = 0) 

= p*Sa (^1 + — j + p*eW,, - p*s\W,^ - p*5W,a - p*s'W^,6a 

-{KQ),a - {'IpS^alp + \\-^P*ma),T + KQW,a + A{KWQ),a\ + 0(4) , (2.14) 



where p is the shear modulus, and 



(2//sf );^ = (2//S6„),5 + - [-(4//iys;.„),;. + 4piy,,s J + 0(4) . 



(2.15) 



The shear modulus p and the compression modulus k are related with the Lame parameter 
A by the relation 

«; = A + 2p/3. (2.16) 
The shear tensor Sgb and the volume dilatation are given by 



Sab 



1 + 



2W 



kOab 



+ 0(4), 



(2.17) 



e = s\k + 4 {^W,ks'' + 3SW) = V-s + ^{VW ■s + SW) + 0(4) . (2.18) 



The gravitational scalar potential W and the vector potential W"^ can be decomposed into 
a sum of two contributions (see paper [5]) 

W^W++W, Wa = W+ + Wa, (2.19) 

where W'^ and are the self parts (resulting from the gravitational action of the body 

under consideration) and W and Wa are the external parts (describing tidal and inertial 
forces) of the metric potentials. Since we take the elastic mechanical ground state as an 
isolated body the external potentials vanish for this state 

W ^Wa^O, W^W+, Wa = W+. (2.20) 

Similarly for the perturbed state we write 

5W = 6W+ + 6W 5Wa = 6W+ + 5W„ . (2.21) 

The external parts of metric potentials, {SW, 5Wa), that result from the ephemerides of 
external bodies are assumed to be known, so only the Eulerian variations of the internal 
potentials have to be determined self-consistently by partial differential equations. The 
perturbed Euler equation can be rewritten as 



+KeW+ + A{KW+e),a]: (2.22) 

where 

r^SW,a + ^(5Wa),T (2.23) 

is the external tidal force density. 

To facilitate the comparison with well known results from Newtonian theory we write 
the equations in a form that was employed in Ref. [9]. To this end we first introduce the 
post-Newtonian Cauchy elastic stress tensor T"'': 

T"^ = KQ6ab + 2pSab - iyl^Sab + nQSab) + 0(4) . (2.24) 

Then Eq.(2.22) takes a form that can easily be compared with the results from Wahr's paper 
(Ref.[9]) 

1^1 + p*s = -[p*eVW+ - p*VW+ ■ {Vsf - p*VSW^ - p*s ■ V{VW+) 



1 



+v • r + ^ |-«:eviy+ + 2iivw 

+^p*(5W+),T + p*f + 0(4), 



(Vs + (Vs 



-V-s / 



(2.25) 



where / is the second rank identity tensor and the superscript T denotes the transpose (not 
to be confused with the local time variable). 

Let us now introduce the following notation generalizing the one from [9] : 

W'^{r) — — $(r) the potential of the ground state 

5W^ = —(j)f the incremental Eulerian gravitational potential energy 

SW = —(pT the tidal potential 

{SW~^°')^T =A ° the time derivative of the self part of vector potential 

{SW^)^T = the time derivative of the external part of vector potential 

f = W5W + iiSW),T = -V0T + 

the tidal force 

(2.26) 

Going into Fourier space with local time variable T and replacing d/dT by iuo the first 
perturbed field equation (Eq.(2.12)) becomes 

^V^ + ^) = 47rG'{ V ■ (ps) + ^ [V ■ (ps) + 3ps • VW+ + p{6W^ + 5W) 

+2V • (pW+s) + 3kV • s] } + 0(4) , (2.27) 

since the tidal potential 5W satisfies D'Alembert's equation. 

The second perturbed field equation (Eq.(2.13)), for the same reason, takes the form 

V'(5iy+"),T = V' 1 " = 47100; + 0(2) . 
By using the definition above (Eqs.(2.26)), the Eq.(2.27) can be rewritten as 

- (v^ + ^) '/•f = 47rG{ V • (p*s) + ^ [-5p*s • V$ - p*0f - p>T 

-2$V • (p*s) + 3kV • s] } + 0(4) . 

The perturbed Euler equation (Eq.(2.25)) takes the form 
4$\ 



(2.28) 



(2.29) 



1 - ^ p*u;'s = - p*V0f + p*V* • (Vs)' + p*s ■ V(V$) 



(Vs + (Vs)^ 



-V • s / 



+^P*{A +A)-p*V0T + O(4), 



(2.30) 



where 

T"^ = KQSab + 2fXSab + [f^Sab + Ke5ab) + 0(4) , (2.31) 

and 

e = V-s-^(v*-s + 0f + 0t)+O(4). (2.32) 

Eqs.(2.28), (2.29), (2.30) and (2.31) are the ones that will now be expanded in terms of 
generalized spherical harmonics (scalar-, vector- and tensor spherical harmonics). 



III. GENERALIZED SPHERICAL HARMONICS 



Eqs.(2.28) — (2.31) are complicated partial differential equations (PDE) that will be 
turned into a set of ordinary differential equations (ODE) by means of expansions in terms 
of scalar-, vector- and tensor spherical harmonics or, briefly, in terms of generalized spherical 
harmonics that are described e.g., by Phinney & Burridge [14], Smith [15] and Wahr [9]. 
Here tensors of arbitrary rank are represented by their components along the complex basis 
e_, eo, e+, where 

e_ = -^(eg - ie^) , eo = , e+ = --^(e^ + ie^) (3.1) 

and e^, eg and e,^ are (Euclidean) unit vectors in r, 9 and direction. Components of a 
tensor field along the complex basis are called canonical components and the basis itself is 
called canonical basis. Generalized spherical harmonics (GSH) are denoted by D^„(^, 0) 
and defined by 

Dl^ie, 0) = (-l)-+-Pr(cos^) exp(im0) , (3.2) 
where P""* is a generalized associated Legendre function 



-1) 



l—n 



2^(1 -ny. 

l—m 



{l-ny.{l + m)\ 

{i + ny.{i-my. 



;i_2;)("-™)/2(l + 2;)- 



■(m+n)/2 



X 



dx . 



l-x)'-"(l + x)'+" . (3.3) 



L>^^ are defined only for / > and \m\ < I (sec, e.g., [14] and [15] for a complete discussion). 
By using GSH, a scalar, vector and 2nd rank tensor can be expanded as 

oo I 

cj>{r,9,cj>) = E E ^Tir)DUO,^), (3-4) 



+1 oo I 

u(r,^,0) = E E E «rW^L(^,0)en, 

n=— 1 1=0 m=—l 

+1 oo I 

T{r,e,<p) = E E E ^r'( 

a,b=—l 1=0 m=—l 



(3.5) 
(3.6) 



The are proportional to the usual spherical harmonics (-D^q — {~^)^\J^'^/ (2^ + 1)^;"*) 
and the D^rnA^A) are similar to the Y;"*"(^,0) of paper [14] (Dj„„ = (_i)m+ny^mn) g^^^g 
useful formulas for objects involving GSH can be found in the Appendix. 



IV. POST-NEWTONIAN PERTURBATION EQUATIONS 

In this section we will rewrite Eqs.(2.28), (2.29), (2.30) and (2.31) in a compact form by 
using GSH. First we expand the scalars, vectors and tensors appearing in these equations 
in the following way 

oo / 



0f(r,e> 



A\r, 9, 4> 



E E ^T{r)DUOA), 

1=0 m=-l 
oo / 

E E ^mnUoA) 



1=2 m=-l 
oo I 

E E srir)DUO,<P), 

1=0 m=-l 
oo I , 

E E Arir)DLiOA), 

1=0 m=-l 
oo I 

Y.Y.X\r)DUeA). 

1=0 m=-l 
oo I 

E E Tr\r)Di^a^,^{eA). 

1=0 m=-l 



(4.1) 
(4.2) 
(4.3) 
(4.4) 
(4.5) 
(4.6) 



Here s" denotes the displacement field, T"^ the incremental Cauchy elastic stress tensor 
(a, b = —1, 0, +1). The gravitational potential $ of the ground (equilibrium) state should 
not be expanded, since it is a known function of the radial coordinate r. 

A. Field equations 



By using the expansions from Eqs.(4.1) — (4.3), the formulas given in the Appendix and 
the orthogonality properties of for each appropriate value oil, m and n, the Eq.(2.29) 
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reduces to 



-47rG 



+ 



(4.7) 



To write this equation in a more compact form it is useful to define new scalar functions 
([9]) 

— 



wr 



> . 



Let 



then 



pm r-pmOO 

rpmO+ _j_ rpmO- 

pm rpmO+ rpmO- 



(4.8) 



(4.9) 
(4.10) 



Substituting Eqs.(4.8), (4.9) and (4.10) into Eq.(4.7), we get 



d_ 
dr 



9i 



\ c 
+(2p*$ - 3ft;) 



,2 



d 



Li. 



dr r J r 



+2t(lp-)c/r|. 



(4.11) 



where 



90 = 

dr 



1 1(1 + 1) 



(4.12) 
(4.13) 



In the Newtonian limit Eq.(4.9) and Eq.(4.11) are correspond to Eq.(III.33) and Eq.(III.34) 
of [9] respectively. 
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Next we consider the second perturbed field equation Eq.(2.28). Expand A ^ ^'^^'^^ ac- 
cording to Eq.(4.4), i.e., 

oo I 



(4.14) 



=0 m=-l 



For the exphcit calculation we have 



+ 



V-VA =C- 



(4.15) 



+ 



Since V A is not a symmetric second rank tensor so that we cannot use Phinney's formulas 
directly. We obtain (for more details see the Appendix) 

rf2 2d 2 



oo / 

e° - E E 

/=0 m=-l 

oo I 

r-EE 

1=0 m=-l 

oo I 

1=0 m=-l 



dr^ r dr r 



™° -^(5r+3r)U!no, (4.16) 



d^ 2d 2 , „ ■ 
ar^ r dr 



(4.17) 
(4.18) 



where Li is defined in Eq.(4.13), and 



'(/-l)(/ + 2) 



(4.19) 



The expansion of takes the form 



CO I 

: 1 

=0 m= 



(4.20) 



Using this expansion in the field equations (2.28) and the orthogonahty relations of 
one finds 



' d^ 


2 d 

r dr 


dr'^ 


''d^ 


2 d 

r dr 


dr'^ 


''d^ 


2 d 


dr^ 


r dr 



2 



2Li 



r 

4Li 



2 ET-'^^Gco'pUr, 

2 ^Jjtrm 



BT 
ET 



(4.21) 
(4.22) 
(4.23) 
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with the definition 



A r++ A r 



+ 



m+ 



AT 

bT 

—rm- 



--El 
--Hi 



AT - 



it: 



(4.24) 



(The definitions for S^*, -^i" ^e used in Section IV.C.) For numerical apphca- 

tions it is useful to reduce the set of seconder order differential equations (Eq.(4.21) — (4.23)) 
to an equivalent set of first order. With 

d 



ET 



dr 
d 

dr 

the post-Newtonian field equations, Eqs.(4.21) - (4.23), finally take the form 

is ^ + ^(1 + i^?) s r + ^ ^ r + ^^Gu^'pur , 



±eT ^ --eT + -,{1 + Ll) ET +^ ' 
dr r r^^ 



BT + ^T^Guj^pVT" , 



dr^' 



r 



^{1 + Ll) H T + ^^Guj'^pWr . 



(4.25) 
(4.26) 
(4.27) 

(4.28) 
(4.29) 
(4.30) 



Compared with the Newtonian case we have six additional equations (Eqs.(4.25 — 4.30)) and 
six additional unknown functions bT, eT, hT, B Ti E T^ H T ^^^^ o^^Y appear in the 
post-Newtonian formalism. 

B. Stress-strain relation 



According to Eqs.(2.31), (2.32) and (2.17), and the stress tensor takes the form 



r = AV-s I +11 Vs+ (Vs) 



+ 



K 



4$V-s - 3(V$ • s + 0f + 0t) 



Using the GSH expansion for the Cauchy elastic stress-tensor one finds that 

oo I 



1=0 m=-l 



= AV-se"^ + pis"^^ + «''■") + 4 4$V-s - 3(V$ • s + 0f + 0t) e"^ , (4.32) 



(4.31) 



11 



where Cab and e"^ are the canonical components of the identity tensor and are equal to 



( 0-1^ 

1 

V-1 y 



a, 6 = -, 0, + 



Using 



oo I 
i=0 m=-l 



ar r r 



1 + 



D 



mO 



and the orthogonality properties of D^^^ we get 



rpm 



rpmO— 

rpmOO 



r 



mO 



mO 



3k / (i 
\ dr 



rpm-\ — 



^^m — I 



rpmO+ rpm+0 ^ 

rpm+ + 



(ir r 



L 



1 ^mO 

"5; 



r 



The relation for T["^ = Pf^ can be written in the form 



(A + 2/x) + — $ 



(4.33) 



(4.34) 
(4.35) 

(4.36) 
(4.37) 



(4.38) 



(4.39) 
(4.40) 
(4.41) 



(4.42) 



which is the post-Newtonian version of Eq.(II1.35) from Wahr (1982). 
By means of Eqs.(4.37) and (4.40) we get 

fi^vr ^QT+ ^(vr + 2L,un , (4.43) 

n^wr = RT + -wr , (4.44) 

dr r 

that generahzes Wahr's Eqs.(III.36) and (III.37) to the post-Newtonian level. Eqs.(4.42), 
(4.43) and (4.44) are the post-Newtonian version of the stress-strain relation. 
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C. Euler equation (Conservation of momentum) 



According to Eqs.(2.30) and (2.32), in the complex basis e_, eo, e+ the perturbed Euler 
equation can be expressed as 



(V-s) - - {^''s^e,a + 0f + 0r)l [ + - pVt" + -^P*{A " + ^") 



+- «;$'"(V-s)-2/x 



(4.45) 



where C" = (v- = T'^^'^ebc- This relation comprises the three equations for a — 
(- 0,+). 

By means of Eqs.(4.36) and (4.39), the a — — equation can be written in the form 



^ rpmO- 



3. 
r 



Li j fi + X 
r I r 



Li(s™+ + s™-) + 2s™° 



r 



I r 



3k 
r 



dr r ) ^ 



r 



(4.46) 



The a = equation, by using Eq.(4.39), can be written as 



' r ' r ^ ' ^ ^ r] r 



+4< 



— $ - (A + 2/^)^0 



r 



4 
r 



p*^o^r - ^P*9o{sr + sT-) + P^ r^r) 



— + - s ' + 
ar r I 



— $ - (A + 2/i)5o 
r 

2k, 



+4$p*a;2sl"° + 87rG'(p - p*$) ^ , 



(4.47) 
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where we have used Eq.(2.6). The pressure p is determined from the equation of state 
P — P{p) that is assumed to be a known function inside the reference body. 



Finally, the a — -\- equation can be brought into the form 



^ 7->mO+ 



L,{sT^ + sT-) + 2s, 



mO 



r 

+ -<! — $Li 



r 



3k 
r 



dr r r 



vu r+A ) ■ 



(4.48) 



We will also rewrite these PN Euler equations to facilitate the comparison with the 
corresponding Newtonian equations. With the definitions from Eqs.(4.8), (4.9) and (4.24), 
the a — Equation (4.47) can be rewritten as 



^ pm 

dr^ 



IT 2 

—pr - -QT + - 

ip rjrt rjrt 



r ^ ' ^ ' dr ^ 



- (A + 2ii)go 

2k\ 



dr r r 



d 



+3 {p*g^ - ^) {g^Ur + C + ""O + 4$p*u;^t/r + 87rG(p - p*$) 
-Ap\BT+X)], (4.49) 



generahzing Wahr's Equation (III. 30). 

Similarly, combining Eq.(4.48) with (4.46) we get 



d_ 
dr 



QT 



3 

--( 

r 



2Li 



2Li 



r 



1 I S'^^r 



^Li(^oc/r+c+^c) 



2Li 



+ A^p'-u'Vr - Ap\E T + K) , (4.50) 
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which is the post-Newtonian version of Wahr's Eq.(III.31). 

Wahr's Eq.(III.32) is generahzed if we combine Eq.(4.48) with (4.46) (using Eq.(4.44)) 

±Rr ^-il-f) RT + '-^wr - (i - ^) ,.^wr - ^(5 r + ^) • 

(4.51) 

Up to now, we have fourteen equations (Eqs.(4.9), (4.11), (4.25), (4.26), (4.27), (4.28), 
(4.29), (4.30), (4.42), (4.43), (4.44), (4.49), (4.50) and (4.51)) for the determination of 
fourteen unknown functions (C/^, F^, W^, g^^, (l)^,, P}^, 

Qmi R-mi ) 

in contrast to eight equations and eight unknown functions in the Newtonian hmit (Ref.[9]). 
V. CONCLUSION 

In this paper we have successfully derived the post-Newtonian equations for the dynamics 
of some elastic deformable astronomical body expanded in terms of Generalized Spherical 
Harmonics(GSH). We have shown that in the Newtonian hmit they reduce to the corre- 
sponding Newtonian equations that can be found in Wahr's paper[9]. The importance to 
separate a set of partial differential equation (PDE) into a set of ordinary differential equa- 
tion (ODE) is obvious for numerical calculations. In that form the equations of motion can 
be directly used in numerical computer programs so that the orders of magnitude of the 
various relativistic terms for concrete problems can be assessed. 

In this paper we considered a non-rotating body as reference state for our perturbation 
theory. Since the precessional and nutational motion of some astronomical body depends 
upon its dynamical ellipticity that is related with its rotational motion due to the inertial 
forces, this formalism cannot be applied directly to problems of precession and nutation. 
The derivation of corresponding equations and junction conditions represented in terms of 
GSH for a rotating elastomechanical ground state will be the subject of a forthcoming paper. 

Acknowledgments 

This work was supported by the National Natural Science Foundation of China (Grant 
No. 10273008) and the German Science Foundation (DFG). 



15 



APPENDIX: SOME USEFUL FORMULAE FOR TENSOR OPERATIONS 



Based on the results of Phinney's paper [14], we derive here some useful formulas for 
computing the perturbed PN equations. The expansion of scalars, vectors and tensors is 
given in Eqs. (4.1)-(4.6). For simplicity, in the following we will omit the summation symbols 

;=0 2^m=-l- 



1. Gradient and Laplacian of a scalar ^{r,0,t 



By means of GSH, a scalar can be expanded as 
and the gradient is given by 



where 



The Laplacian of $ can be written as 



(f ^2 d + 1) 



(A.l) 
(A.2) 

(A.3) 



dr"^ r dr 

where Ca/s are the canonical components of the identity tensor, its value is given by Eq.(4.33). 
For simplicity, we represent [V(V$)]"'^ as <l>'"^, then the second partial derivatives of $ are 
given by 

1 



I 



r dr 



mO 



m+2 ) 



(A.4) 
(A.5) 

(A.6) 

(A.7) 
(A.8) 
(A.9) 
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where 



2. Operations involving a scalar $(r) 

If a scalar is a function of the radial coordinate r only, i.e., $ = ^{r), we have 

= 0, 

V2$- (— + 1—] $ 
r or 

$,00 ^ 

3. Gradient and divergence of a vector-field 

The expansion of a vector is given by Eq.(4.3), i.e., 

u-{r, e, 0) = Ur{r)Dl,{e , 0) . (a = -1 , , +1) 

We write the vector gradient as 

where 



' -i [L^ur + e«_(-i)"t/r° + Soai-irur) , = -i 
TrUr, p^o 



where 



Lq — Li , L_i — L2 
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The symmetric vector gradient can be written as 

-1 r. 



2r 



d: 



mO ' 



_^ 1 I jjma _ -^1 jjmO 

,dr r / ' r ^ 



2 ^ ^2 
The divergence of a vector-fields is expanded in the form 

d 2 



{a = ±1) 



(A.22) 
(A.23) 



V • u = u'^'f^ea 



dr r 



mO 



1 



D 



mO 



(A.24) 



4. Divergence of a symmetric, second order tensor-field 

For a symmetric tensor-field T"'^, the divergence, (V- T)" = T"'^''''e^^ = C", is given by 
2^ 



C° = 



" , - 1 T^mOO , , 1/^(^ + 1) 



mO+ _|_ rprnO 



d 3' 

dr r 



(A.25) 



■ (« = ±1) 



(A.26) 



5. Laplacian of a vector-field 

The Laplacian of a vector-field can be written in the form 

(V^A)° = (A^'^^'T + A'^'"'^ - A^'^'°)e;3^ = r , (a = -l, 0, +1) (A.27) 

where A""'^ is an unsymmetric second order tensor. In terms of GSH, a tensor-field of second 
rank can be expanded as 

m-^ir, e, 0) = Mr^r)Dl,^,^f,)(0, 0) . (a, /3 = -1, 0, +1) 

Differentiation yields 

where the calculation of the quantity M'^"'^^'^ is similar to Eq.(2.13) of Ref.[14], note that 



Di 



Some terms, used in this paper, read 



mOO|0 



dr ' 



(A.28) 
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mO-|+ 



mOO^ 



Ml 



mO+\- 



M, 



TO-0|0 



Ml 



m 1 + 



M, 



m-+\- 



m+0|0 



Ml 



m+-\ + 



M, 



m++\- 



Now we have 



Set 



then 



m-f'ir, e, 0) = A"''^ = Mr''(r)i?^(„+^)(^, , (a, /3 = -1, 0, +1) 



(A.29) 
(A.30) 
(A.31) 
(A.32) 
(A.33) 
(A.34) 
(A.35) 
(A.36) 



(A.37) 



The quantities M;"'"^l^(r) are derived from Eqs.(A.28), (A.29), (A.30), (A.31), (A.32), 
(A.33), (A.34), (A.35) and (A.36), and M;"*"^(r) = ^^'"'^(r) can be calculated by Eq.(A.20), 
then we can deduce 



r 





2 d 




r dr 


''d^ 


^2 d 
r dr 


dr"^ 


''d^ 


^2 d 
r dr 


dr^ 



2Li 



A] 



r 

- 4'' 



(Ar + Ar) Di, , 



mO (. T^l 



Ar^ 



2Li 



m+1 ■ 



(A.38) 
(A.39) 
(A.40) 
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